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Abstract 

Background: Action potentials can be initiated at various subcellular compartments, such as axonal hillock, soma 
and dendrite. Mechanisms and physiological impacts for this relocation remain elusive, which may rely on input 
signal patterns and intrinsic properties in these subcellular compartments. We examined this hypothesis at the 
soma and axon of cortical pyramidal neurons by analyzing their spike capability and voltage-gated sodium channel 
dynamics in response to different input signals. 

Results: Electrophysiological recordings were simultaneously conducted at the somata and axons of identical 
pyramidal neurons in the cortical slices. The somata dominantly produced sequential spikes in response to long-time 
steady depolarization pulse, and the axons produced more spikes in response to fluctuated pulse. Compared with the 
axons, the somata possessed lower spike threshold and shorter refractory periods in response to long-time steady 
depolarization, and somatic voltage-gated sodium channels demonstrated less inactivation and easier reactivation in 
response to steady depolarization. Based on local VGSC dynamics, computational simulated spike initiation locations 
were consistent with those from the experiments. In terms of physiological impact, this input-dependent plasticity of 
spike initiation location made neuronal encoding to be efficient. 

Conclusions: Long-time steady depolarization primarily induces somatic spikes and short-time pulses induce axonal 
spikes. The input signal patterns influence spike initiations at the axon or soma of cortical pyramidal neurons through 
modulating local voltage-gated sodium channel dynamics. 
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Introduction 

Neurons integrate synaptic inputs and produce sequen- 
tial spikes as digital codes in the brain. In terms of the 
sources of producing spikes, previous studies indicated 
that a single spike was initiated at axonal hillock [1-10], 
as well as somata and dendrites [11-21]. These data raise 
the question whether the location of spike initiation 
undergoes plasticity, in which input signals and local in- 
trinsic properties may be involved. The elucidation to 
the dynamical alternation of spike initiation locations is 
critically important to understand how the neurons 
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integrate synaptic signals and produce their digital codes 
in the brain efficiently. 

Synaptic input signals that evoked neuronal spikes 
in vivo appeared long duration [15,22-25]. This long- 
time physiological depolarization induced sequential 
spikes dominantly at the somata [15], whereas short 
pulses initiated individual spikes at the axonal hillock 
[2,3,7,9]. Why the long-time versus short-time input sig- 
nals influence the locations of spike initiation remains 
elusive. In addition, long-time signals integrated from 
synaptic inputs in vivo are classified into steady 
depolarization and fluctuated pulses [15]. It needs to be 
addressed whether these two patterns of input signals 
evoke sequential spikes at different subcellular compart- 
ments, an input-dependent plasticity of spike initiation 
location. 
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In terms of mechanisms underlying the subcellular 
localization of spike initiation, the dynamics and density 
of local voltage-gated sodium channel (VGSC) are pre- 
sumably involved. It was suggested that a high density of 
VGSCs is critical for spike initiation at the axonal hillock 
[2,7,9]. However, the depolarization did not increase the 
number of functional VGSCs at axonal initial segment 
[26]. The intact initial segment with dense VGSC clus- 
ters was unnecessary for inducing spikes in the neurons 
[27]. VGSCs dynamics also played an important role in 
spike initiation location [3,28]. Therefore, the changes in 
local VGSC dynamics and/or density are likely involved 
in the plasticity of spike initiation location. As the activity- 
dependent redistribution of high dense VGSCs occurred 
within axonal hillock and took a long time [29,30], the 
plasticity of spike initiation location may be based on local 
VGSC dynamics. In order to address the roles of input sig- 
nal patterns and local VGSC dynamics in spike initiation 
relocation, we analyzed the sequential spikes and VGSCs 
dynamics at the axonal bleb and soma of identical pyram- 
idal neurons simultaneously in sensory cortical slices. 

Results 

Membrane depolarization signals in vivo are long time, 
whose patterns are generally classified into steady and 
fluctuated pulses [15]. The action potentials can be evoked 
at various subcellular compartments, such as axonal 
hillock, soma and/or dendrite [2,3,11-13,15]. We have 
proposed to examine whether input signal patterns 
influenced spike-initiation location as well as how local 
VGSC dynamics regulated this input-dependent reloca- 
tion of spike-initiation sites. In terms of the strategies 
to address these issues, we mainly analyzed correlations 
between input-signal patterns and spike-initiation loca- 
tions by changing input signal patterns and subcellular 
compartment functions. This analysis would indicate 
whether the dynamic relocation of the spike-initiation 
sites was naturally present. To its underlying mecha- 
nisms, we focused on analyzing VGSC dynamics at 
these subcellular compartments in response to different 
input signals. Based on the dynamics of local VGSCs, 
we conducted computation simulation to test whether 
their dynamics characteristics were responsible for spike- 
initiation relocation. Furthermore, we examined whether 
the changes of local VGSC s dynamics would shift spike- 
initiation location. Finally, we studied physiological 
impacts for the spike-initiation relocation, such as the 
efficacy of neuronal encoding. 

Input signal patterns influence spike-initiation at the 
somata and axons 

If steady depolarization versus fluctuated one initiate 
sequential spikes at the different locations of given neurons, 
altering input signal patterns should drive spike-initiation 



relocation. In other words, if one of the subcellular com- 
partments is preferentially sensitive to an input pattern for 
firing spikes, the spike capability in response to this input 
signal should be higher in this compartment than others. 
While changing input signal patterns, we analyzed the 
input-output curves of these compartments to assess 
their sensitivity to input signals and their spike capability. 

Whole-cell recording was simultaneously conducted 
on the soma and axonal bleb of the same cortical pyr- 
amidal neurons to apply different input signals locally 
and to acquire sequential spikes (Figure 1A). Long-time 
(1 second) steady depolarization pulses or fluctuated 
pulses (cosine-wave, [31]) were injected into these two 
compartments to induce spikes, respectively (Additional 
file 1: Figure SI). While these depolarization pulses were 
injected into either of two compartments, the same 
amount of spikes was recorded at these two compart- 
ments. This spiking fidelity was due to the secure spike 
propagation on cell membrane between the axonal blebs 
and somata [13]. On the other hand, the subthreshold 
potential was electrotonically propagated and shunted by 
potassium channels, such that they would be decayed. We 
evaluated the capability of firing spikes (input-output and 
spike threshold) in these two compartments by injecting 
pulses and evoking spikes locally. 

No matter what the axon or soma is a location to initiate 
spikes, this compartment should have higher capability 
to convert input signals into spikes, i.e., more efficient 
input-outputs transformation. Steady pulses in various 
intensities were injected into the soma and axonal bleb 
(20 ~ 50 (im away from the soma), respectively, to induce 
spikes (red trace in Figure IB for the soma and blue for 
the axon). Figure 1C shows spikes per second versus 
stimulus intensities at the soma (red symbols) and axon 
(blue) in this example. Somatic input-output curve (red 
triangles in Figure ID, n = 20) appears on the left-top side 
of axonal one (blue circles, n = 20; asterisks, p < 0.05, 
paired-t test), indicating that somatic spike threshold is 
lower and identical stimuli induce more somatic spikes. 
The somata are sensitive to long-time steady input signals. 

On the other hand, the fluctuated pulses (black trace 
in Figure IE) in various intensities were injected into the 
soma and axon of these pyramidal neurons, respect- 
ively, to induce spikes (red trace for the soma and blue 
for the axon). Axonal input-output curve (blue sym- 
bols in Figure IF ~ G) is on left-top side of somatic one 
(red, n = 20; asterisks, p < 0.05, paired-t test), indicating 
that the axons are more sensitive to the fluctuated in- 
put signals than the somata. 

We examined these indications under a voltage-clamp 
to reduce the effect of passive membrane property on the 
recording. Steady depolarization induced more spikes at 
the soma (red trace in Figure 2 A) than at the axon (blue). 
Somatic input-output curve (red triangles in Figure 2B) is 
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Figure 1 Long-time step pulses initiate sequential spikes at the soma of cortical pyramidal neurons, but the fluctuated signals induce 
spikes at the axon. A) Top panel shows the images of dual recording on the soma and axonal bleb under fluorescent (left)/DIC (right) microscope. 
Bottom panel shows an image of neurobiotin-labeled pyramidal neuron. B) Long-time step depolarization (middle black trace) induces more spikes at 
the soma (red trace) than at the axon (blue). C) shows the number of spikes vs. stimulus intensities at the soma (red symbols) and axon (blues) of this 
pyramidal neuron. D) illustrates input-output curves at the soma (red symbols) and axon (blues; n = 20, p < 0.05). E) The fluctuated signal (a cosine 
wave, middle black trace) induces more spikes at the axon (blue trace) than at soma (red). F) shows the ratio of spikes to cosine-waves vs. stimulus 
intensities at the axon (blue symbols) and the soma (reds) of this pyramidal cell. G) shows input-output curves for the axon (blue symbols) and the 
soma (reds; n = 20, p < 0.05). 



on the left-top side of axonal one (blue circles, n = 9; 
asterisk, p < 0.05, paired-t test). On the other hand, the 
fluctuated signals induced more axonal spikes (blue 
trace in Figure 2C) than somatic spikes (red). Axonal 
input-output curve (blue circles in 2D) appears on left- 
top side of somatic one (red triangles, n = 9; asterisks, 
p < 0.05, paired-t test). The data are consistent with 
those under current-clamp (Figure 1). 

If the soma is sensitive to steady pulse and the axonal 
segments are sensitive to fluctuated pulse, each of these 
compartments would fire the spikes in the highest fre- 
quency when receiving its sensitive signal. This assump- 
tion is based on a principle that sinoatrial node with the 
highest pacemaker spikes controls the heart rate. Steady 
depolarization pulses in various intensities were injected 
into the soma and different axonal segments, respect- 
ively, to induce spikes (Figure 3A). Figure 3B illustrates 
input-output curves for the somata (red symbols, n = 49) 
and axonal segments away from the soma in 5-29 um 
(light-blue, n = 36), 30 ~ 55 um (purple, n = 26) and >55 um 
(dark-blue, n = 7). The maximal number of spikes from 3B 
versus the distances between axonal segments to soma 
were plotted in Figure 3C. The sensitivity to steady 
depolarization and the capability to produce its induced 
spikes are decreasing from the soma toward distal axons. 

Moreover, the fluctuated pulses at various intensities were 
injected into the soma and different axonal segments, 



respectively, to induce spikes (Figure 3D). Figure 3E 
illustrates input-output curves at the soma (red symbols, 
n = 19) and the axonal segments away from the soma in 
5 ~ 29 um (light-blue, n = 13), 30 ~ 55 um (purple, n = 10) 
and >55 um (dark-blue, n = 10). The efficiency of convert- 
ing the fluctuated signals into the spikes (the ratio of 
spikes to cosine-waves; dark-red) versus the distances of 
axonal segments to the soma were plotted in Figure 3F. 
The sensitivity to fluctuated signals and the capability to 
fire its induced spikes are increasing from the soma 
toward distal axons. 

We analyzed input- output curve between the soma 
and axon segments of identical pyramidal neurons by alter- 
ing input signal patterns. Long-time steady depolarization 
preferentially induces sequential spikes at the soma, but 
fluctuated one (sequential short-time pulses) initiates spikes 
at the axons. Moreover, we confirm this input-dependent 
plasticity of spike initiation location by changing com- 
partment functions to shift spike initiation locations. It 
is noteworthy that spike initiation relocation is not 
caused by passive membrane property since there is no 
difference in the input resistance between the soma and 
axon (Additional file 2: Figure S2). As spike initiation is 
presumably controlled by VGSCs, we propose that the 
sensitivities of the somata and axons to the different 
input signals are due to the difference of their local 
VGSC s dynamics, which we examined below. 
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Figure 2 Long-time step pulses preferentially initiate spikes at the soma of cortical pyramidal neurons under the condition of 
voltage-clamp recording, but the fluctuated signals induce spikes at the axon. A) A long-time depolarization voltage (50 mV, middle black 
trace) induces more spikes at the soma (red trace) than at the axon (blue). B) illustrates input-output curves at the soma (red symbols) and the 
axon (blues; n = 9, p < 0.01). C) A fluctuated voltage signal (cosine-wave, 50 mV; middle black trace) induces more spikes at the axon (blue trace) 
than at the soma (red). D) shows input-output curves for the axon (blue symbols) and the soma (reds; n = 9, p < 0.01). 



Somatic and axonal VGSCs are different in response to 
long-time and short-time signals 

The VGSC dynamics sets spike capability as well as thresh- 
old potentials and refractory periods [32,33]. In examining 
the role of somatic and axonal VGSCs in this input- 
dependent plasticity of spike initiation location, we ana- 
lyzed the influences of input signal patterns on spike 
thresholds and refractory periods at the soma versus the 
axon of pyramidal neurons. We also examined the influ- 
ences of changing VGSCs dynamics on spike thresholds 
and refractory periods at the soma and axon. 

Somatic and axonal spike thresholds and refractory periods 
change in response to different signals 

Spike thresholds and refractory periods were measured 
at the soma and axonal bleb of identical cells [13]. In the 
measurements of somatic and axonal thresholds, pulse 
durations and intensities were inversely adjusted to 



induce the spikes at threshold level. Axonal spike 
thresholds appear low in response to short-time pulses 
(blue trace in Figure 4A), but somatic thresholds are 
low to long-time pulses (red). Threshold stimuli versus 
depolarization time at the soma (triangles) and the 
axon (circles) are plotted in Figure 4B, in which their 
fitting curves cross at 20 ms (asterisks, p < 0.05, n = 12). 
In addition, the refractory periods of axonal spikes 
induced by short-time pulses appear short (blue traces 
in Figure 4C), while those of somatic spikes by long-time 
pulses are short (reds). Figure 4D shows spike refractory 
periods versus depolarization time at the axon (circles) 
and soma (triangles, n = 9; p < 0.05). When long-time 
depolarization signals are inputted to pyramidal neu- 
rons, somatic spike thresholds and refractory periods 
are converted to be lower than axonal ones, so that 
sequential spikes are initiated primarily at cell body, or 
vice versa. 
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Figure 3 The axon possesses the maximal ability of firing spikes in response to fluctuated signals, and the soma has the maximal 
ability of firing spikes in response to long-time steady depolarization. A) An example shows spike waveforms from the soma to different 
axonal segments. B) shows the input-output curves for the soma (red symbols) and different axonal fragments in 5 ~ 29 urn (light-reds), 30 ~ 55 urn 
(pinks) and >55 urn (blues) away from the soma, in which sequential spikes are induced by long-time steady depolarization. C) illustrates the maximal 
values of input-output curves (normalized stimuli at 10) for steady-state pulses versus distances of axonal bleb to the soma. D) A sample shows spike 
waveforms from the soma to different axonal segments. E) shows input-output curves for the soma (red symbols) and distinct axonal fragments in 
5 ~ 29 urn (light-reds), 30 ~ 55 urn (pinks) and >55 urn (blues) away from the soma, in which the spikes are induced by the fluctuated depolarization 
(cosine wave at 50 Hz). F) shows the maximal values of input-output curves for the fluctuated pulses vs. distances of axonal bleb to the soma. 



We also examined the effects of reducing VGSC in- 
activation on spike threshold and refractory period. The 
reduction of VGSC inactivation was fulfilled by using 
anemone toxin (ATX), a blocker of VGSC inactivation 
[34,35]. ATX appears to reduce axonal spike threshold 
(green traces in Figure 5A). Figure 5B illustrates thresh- 
old stimulations versus depolarization time at the axon, in 
which the spike thresholds are attenuated by ATX (green 
symbols; p < 0.05, n = 14). Furthermore, ATX appears 
to shorten axonal refractory periods (green trace in 
Figure 5C). Figure 5D shows spike refractory periods versus 
depolarization time at the axon, in which ATX shortens 
spike refractory periods mainly induced by long-time pulses 
(green symbols; p < 0.05; n = 13). It is noteworthy that ATX 
does not influence somatic spike threshold and refractory 
period (Additional file 3: Figure S3). These data indicate 
that long-time depolarization signals mainly inactivate 
axonal VGSCs and lower its spike capability. 

In terms of the influence of input signals on VGSCs 
dynamics, long-time steady depolarization mainly inacti- 
vates axonal VGSCs, such that the steady signals initiate 



sequential spikes at the somata and the fluctuated signals 
to initiate spikes at the axons. We further investigated how 
long-time steady depolarization pulses influence the inacti- 
vation and reactivation of axonal VGSCs vs. somatic ones. 

Axonal VGSCs are easily inactivated and less reactivated by 
long-time steady depolarization 

VGSCs were inactivated in voltage-/state-dependent 
manners [36,37]. In vivo signals induced sequential 
spikes on the long-time depolarization [13,22-25]. This 
pre-depolarization may inactivate axonal VGSCs [38] to 
elevate spike threshold and refractory period (Figure 4), 
such that spike initiation shifts to the soma [15,39]. To 
examine this implication, we compared VGSC inactivation 
and reactivation at the soma versus axon by cell-attached 
recording simultaneously at the two subcellular compart- 
ments of identical pyramidal neurons (Figure 6). 

VGSCs' inactivation curves were measured by giving pre- 
depolarization pulses ahead of a fixed depolarization pulse. 
Compared to somatic VGSCs (red traces in Figure 6A), 
axonal ones (blue trace) appear obvious inactivation in 
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Figure 4 Spike thresholds and refractory periods are changed dynamically by the patterns of input signals at the soma vs. axon. A) 
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stimuli vs. depolarization time at the soma (triangle symbols) and the axon (circles, n = 12). C) Left panels show that spike refractory periods ap- 
pear short at the axon (blue traces at top panels) by short-time pulses. Right panels show that refractory periods are short at the soma (red lines 
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Asterisks present p < 0.05. 



response to these pre-depolarization pulses. Axonal 
VGSC inactivation curve (blue line/circles in Figure 6B) 
is on the left side of somatic one (red triangles), that is, 
the pre-depolarization pulses make axonal VGSCs to be 
easily inactivated. This result is consistent to a report 
that axonal Navl.6 inactivation is more than somatic 
Navl.2 inactivation during pre-depolarization [40]. 

To test whether the inactivation of axonal or somatic 
VGSCs affects their subsequent activation, we measured 
their activation curves under a pre-depolarization. 
Figure 6C ~ D shows somatic VGSCs activation with 
(top panel in 6C and filled triangles in 6D) and without a 
pre-depolarization (control; bottom in 6C and open trian- 
gles in 6D). Figure 6E ~ F shows axonal VGSCs activation 
under partial inactivation (top panel in 6E and filled circles 
in 6F) and control (bottom in 6E and open circles in 6F). 
Axonal VGSCs become difficultly reactivated if a pre- 
depolarization induces their partial inactivation. 



In summary, long-time signals dominantly suppress 
the reactivation of axonal VGSCs and lead to the high 
values of axonal spike thresholds and refractory periods 
(Figures 4, 5 and 6). So, long-time steady signals mainly 
initiate sequential spikes at the soma, and the fluctuated 
ones induce spikes at the axon (Figures 1, 2 and 3). We 
examined this indication by changing axonal VGSC dy- 
namics and computational simulation. 

The role of local VGSCs in input-dependent relocation of 
spike initiation 

The data above indicate that spike onset at the axon or 
the soma is related to their local VGSC dynamics. To 
examine the role of local VGSCs in the input-dependent 
plasticity of spike initiation location, we conducted 
computational simulation as well as changed local 
VGSC dynamics to see spike initiation relocation. 
Introducing VGSCs, which possess somatic featured 
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VGSC dynamics, into computational modeling is ex- 
pected to simulate sequential spikes similar to somatic 
spikes, or vice versa. The experimental upregulation 
and downregulation of local VGSC dynamics are ex- 
pected to cause the relocation of spike initiation. 

Easier inactivation and less reactivation of axonal VGSCs 
simulate spikes being somatic origin 

In computational simulation, we introduced the curves 
of axonal VGSC reactivation (blue trace in Figure 7 A) 
and somatic one (red) under pre-depolarization (from 
Figure 6D and 6F) into NEURON model [41,42]. In the 
modeling, threshold potentials (AV) for axonal VGSCs 
were more positive than somatic VGSCs, and axonal 
VGSC reactivation (Ai) was 50% lower than somatic one. 
By introducing these factors, long-time steady signal 
induced more simulated spikes at the soma (red trace 
in Figure 7B) than the axon (blue). The threshold of 



sequential spikes is lower, and the number of spikes by 
identical stimuli is higher at the soma (red symbols in 
Figure 7C) than the axon (blue). 

On the other hand, we inputted the curves of axonal 
VGSC activation (blue curve in Figure 7D) and somatic 
one (red) under no pre-depolarization (control in Figure 6D 
and F) into NEURON model. In the modeling, threshold 
potential (AV) and VGSC activation (Ai) for the axon were 
similar to those for the soma. The fluctuated signals under 
this condition induced more simulated spikes at the axon 
(blue trace in Figure 7E) than the soma (red trace). The 
thresholds of sequential spikes are lower, and the number 
of spikes by identical stimulus is higher at the axon (blue 
symbols in Figure 7F) than the soma (red). 

The data from computational simulation supports the 
experimental data that long-time steady depolarization 
mainly inactivates axonal VGSCs and initiates sequential 
spikes at the soma, whereas fluctuated signals initiate 
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Figure 6 Pre-depolarization inactivates voltage-gated sodium channels (VGSC) partially and weakens their subsequent activation at 
the axon. A-B) Pre-depolarization pulses inactivate VGSCs mainly at the axons (A). The inactivation curve of axonal VGSCs (blue symbols in B) is 
on the left side of that of somatic ones (red triangles). Axonal VGSCs are more easily inactivated. C-D) A pre-depolarization pulse weakens somatic 
VGSCs activation slightly (top panel in C for waveform and red-filled symbols in D) compared with control (bottom panel in C and opens in D). 
E-F) A pre-depolarization pulse significantly weakens axonal VGSCs activation (top panel in E for waveform and blue-filled symbols in F) compared 
with control (bottom in E and opens in F). 



the spikes at the axon, i.e., input-dependent plasticity of 
spike initiation location. 

The manipulation of axonal VGSCs function changes the 
location of spike initiation 

Long-time steady depolarization initiates more somatic 
spikes than axonal spikes (Figures 1, 2 and 3) by inactivat- 
ing axonal VGSCs (Figures 4, 5, 6). If it is right, axonal 
spike capability in response to steady depolarization 
should be enhanced by upregulating axonal VGSC dy- 
namics. While puffed ATX to the axonal blebs of pyr- 
amidal neurons (a green electrode in Figure 8A), we 
observed that ATX increased the number of spikes in- 
duced by steady pulses at the axon (Figure 8B). Figure 8C 
shows input-output curves for the soma (red symbols) 
as well as axon before (blues) and after using ATX (greens; 
p < 0.01, n = 10). The prevention of axonal VGSC inactiva- 
tion strengthens axonal spike capability. 

On the other hand, we inactivated axonal VGSCs by 
long-time depolarization pulse to examine whether axonal 
spike capability in response to fluctuated signal would be 
weakened. In dual recording on the soma and axonal bleb 
of identical pyramidal cells (Figure 9A), we inactivated 



axonal VGSCs by using long-time steady depolarization 
pulses in two methods. The fluctuated signals were 
injected into the axons to evoke the spikes, while the 
steady signals were injected into the somata to indir- 
ectly inactivate axonal VGSCs (Figure 9B). The fluctu- 
ated signals were injected into the somata to evoke the 
spikes, while the steady signals were injected into the 
axons to directly inactivate axonal VGSCs (Figure 9C). 
Both approaches reduce spike capability induced by the 
fluctuated signal, but the direct inactivation of axonal 
VGSCs more dominantly reduces spike capability (gray 
symbols in Figure 9D) than the indirect inactivation of 
axonal VGSCs does (orange symbols, n = 7). Neuronal 
spike capability and patterns by inactivating axonal VGSCs 
are similar to somatic spikes (please compare to Figure 3). 
Thus, an inactivation of axonal VGSCs weakens axonal 
spike capability in response to the fluctuated signals, or 
shifts spike initiation location toward the soma. 

The input-dependent plasticity of spike initiation location 
increases neuronal encoding efficiency 

In terms of physiological significances for the plasticity of 
spike initiation location, we assume that this relocation of 
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Figure 7 Computational simulation favors a somatic origin of sequential spikes under the condition of pre-depolarization, and axonal 
origin for spikes without pre-depolarization. A) A pre-depolarization pulse significantly inactivates axonal VGSCs, including an increase of their 
threshold by AV and a decrease of activation by Ai. B) The number of simulated spikes appears higher at the soma (red trace) than the axon (blue) 
induced by a long-time step depolarization (black). C) illustrates spike number vs. normalized stimuli at the soma (red triangles) and the axon (blue 
circles; n = 20, p < 0.01). D) shows the activation curves of somatic VGSCs (red trace) and of axonal ones (blue) by the pulses without pre- depolarization 
and hyperpolarization. E) The number of simulated spikes appears higher at the axon (blue trace) than the soma (blue) induced by a cosine-wave 
depolarization pulse (black). F) illustrates input-output curves at axonal spikes (blue circles) and somatic ones (red triangles; n = 20, p < 0.01). 




Figure 8 ATX dominantly secures the axonal capability of firing sequential spikes induced by long-time steady depolarization pulses. 

A) shows a diagram for whole-cell recording at the soma and axon, and puffing 5 uM ATX (green pipette) toward axonal bleb. B) Long-time steady pulse 
initiates sequential spikes at the axon of cortical pyramidal neurons before (blue trace) and after A~D< application (green), as well as sequential spikes at 
the soma of pyramidal neurons (red trace). C) shows spikes per second versus stimulus intensities at the axons of pyramidal neuron before (blue circles) 
and after ATX applications (greens, n = 10; asterisks, p < 0.05), as well as input-output curves at the somatic recordings (red triangles; n = 10). 
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Figure 9 The reduction of axonal VGSC function by steady depolarization signal may weaken axonal capability to produce spikes in 
response to fluctuated signals. Axonal VGSCs are inactivated by long-time steady depolarization pulses in two ways. A) shows dual recording 
at the soma and axonal bleb. B) The fluctuated pulse (FP) is injected into the axon to induce spikes, while the steady pulse (SP) is injected into 
the soma to indirectly inactivate axon VGSCs. C) FP is injected into the soma to induce the spikes, while SP is injected into the axon to directly 
inactivate axonal VGSCs. D) shows input-output curves for directly (open symbols) and indirectly (filled symbols) inactivating axonal VGSCs. Spike 
capability induced by the fluctuated pulse is more dominantly reduced by the direct inactivation of axonal VGSCs than the indirect inactivation of 
axonal VGSCs does (n = 7). 

V J 



spike initiation may be for the neurons to save the energy 
during neuronal encoding. The neurons have to step over 
two barriers [43], spike thresholds and refractory periods 
[32,44-46], to produce the spikes. If this assumption is 
right, we should see lower energetic barriers for somatic 
spikes by long-time steady signal or axonal spikes by fluc- 
tuated signal. By calculating the energetic barrier, the 
multiplication of spike threshold and refractory period 
[43], we found that this value at the soma was low by 
using long-time pulses (red symbols in Figure 10) and the 
value at the axon was lower by short- time ones (blues). 
Thus, spike initiations at the soma by long-time signals 
and at the axon by short ones make the neurons to utilize 
the lower levels of energy. In other words, the input- 
dependent plasticity of spike initiation location is benefi- 
cial for the neurons to save energy during encoding digital 
signals. 

In addition, the portion of the fluctuated signal was 
about 40% in the integrated synaptic signals [15]. To 
convert synaptic inputs into spikes efficiently, the neu- 
rons should make each of fluctuated signals to induce a 
spike. Where is this optimal conversion fulfilled? As the 
axon is more sensitive to the fluctuated signal (Figures 1, 
2 and 3), it is likely a subcellular compartment that 
converts the fluctuated signals into the spikes with an 
optimal match in their frequencies. We tested this possi- 
bility by measuring the optimal efficiency of converting 
cosine waves into spikes at the soma and the axonal 
segments. 

The fluctuated signals (cosine wave) with frequency 
increments were injected into the soma and different 
axonal segments, respectively, to induce the spikes. If 
the soma or axon in response to these signals reached 



an optimal level (each of pulses triggered a spike), we 
named it as optimal response frequency. The optimal 
frequency of converting fluctuated signals into spikes ap- 
pears the highest at the distal axon (Figure 11 A). Figure 11B 
shows optimal frequency versus distance of axonal seg- 
ments to the soma (n = 17 for soma; n = 43 for axon). 
Figure 11C shows number of spikes versus frequency of 
depolarization (DP) at the axon (blue symbols) and 
soma (reds, n=ll). Therefore, the axonal segments 
away from the soma above 50 \im are the subcellular 
compartment that converts the fluctuated signals into 
the spikes with optimal match in high frequencies, i.e., 
the axons confer the neurons to efficiently encode 
digital spikes in response to the fluctuated input signals. 

Discussion 

Long-time steady depolarization initiates sequential 
spikes dominantly at the somata of cortical pyramidal 
neurons, whereas fluctuated signal instigates the spikes 
at their axons (Figures 1, 2 and 3), i.e., an input- 
dependent plasticity of spike initiation locations. In 
terms of its mechanism, long-time pulses make somatic 
spike thresholds and refractory periods to be lower and 
short-time ones make such parameters to be lower at the 
axon (Figures 4 and 5). Long-time depolarization mainly 
inactivates axonal VGSCs, or vice versa (Figure 6). 
Computational simulation based on the kinetics of axonal 
and somatic VGSCs supports that long-time steady 
depolarization induces somatic spikes and fluctuated 
signals do axonal ones (Figure 7). Our studies in the 
experiments and computational simulation reveal the 
spike initiation relocation between the somata and 
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Figure 11 The axons confer the pyramidal neurons to efficiently encode the digital spikes in response to the fluctuated input signals. 

The fluctuated signals (cosine waves) at threshold stimuli in different frequencies were injected into the soma and different axonal segments, 
respectively, to induce spikes. A) shows the optimal frequency of converting the fluctuated signals into spikes at the soma and axons. B) illustrates 
optimal frequency versus distance of axonal segments to the soma (soma, n = 1 7; axons, n = 43). C) shows the frequency of depolarization (DP) vs. 
number of spikes (n = 1 1). Axonal segments away from the soma larger than 50 urn appear a subcellular compartment of converting the fluctuated 
signals into the spikes with optimal match in their frequencies. 
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axons of cortical pyramidal neurons, which is con- 
trolled by local VGSC kinetics. 

In order to make sure the input-dependent plasticity 
of spike initiation location, we applied multiple strat- 
egies. For instance, subcellular compartments in terms 
of their sensitivity to distinct input signals and their 
ability to produce spikes were analyzed under the con- 
ditions of current-clamp and voltage-clamp record- 
ings. These experiments were conducted from the 
somata to different axonal segments in order to find a 
dominant location of spike initiation. Moreover, the 
methods in the functional upregulation and downreg- 
ulation of subcellular compartments were used to 
analyze the relocation of spike initiation. The results 
proved by these strategies make the conclusion, the 
input-dependent plasticity of spike initiation location, 
to be convincing. 

Current reports indicated that the neuronal activities 
induced the plasticity in terms of the length and distri- 
bution of high dense VGSCs at axonal hillock [29,30]. 
This plasticity of VGSC distribution spends long time 
for its onset and occurs within axonal hillock. Different 
from this slow onset plasticity, the input-dependent 
plasticity of spike initiation location in our study 
develops quickly, shifts between the soma and axon 
(Figures 1, 2 and 3) and depends on local VGSC kinet- 
ics (Figures 4, 5 and 6). This subcellular relocation of 
spike initiation is a novel type of neuronal plasticity, 
compared to long-term plasticity in neuronal excitabil- 
ity [47-54], neuronal homeostasis [55] and VGSC redis- 
tribution [29,30]. 

A major finding in our studies is that long-time steady 
depolarization signals initiate sequential spikes at the 
somata of cortical pyramidal neurons (Figures 1, 2 and 
3) through inactivating axonal VGSCs and conferring 
somatic VGSCs to be less inactivation and easily reacti- 
vation (Figures 4, 5 and 6). As long-time steady pulses 
are a major portion of in vivo signals [15], a somatic ori- 
gin to fire sequential spikes is physiologically important. 
One could argue that the smaller axonal volume than 
somatic one plus the higher density of axonal VGSCs 
may enable the depolarization to bring denser positive 
charges into the axon that facilitates the axons to reach 
a threshold and produce the first spike ahead of the 
soma. However, the voltage-/state-dependent VGSC in- 
activation (Figures 4, 5 and 6; [36,37]) may also make 
this depolarization inactivating axonal VGSCs, espe- 
cially axonal Navl.6 [40], which does not help to fire se- 
quential spikes at the axon. 

On the other hand, the fluctuated depolarization sig- 
nals, a sequence of short-time pulses, initiate axonal 
spikes at cortical pyramidal neurons (Figures 1, 2 and 3), 
due to low threshold and refractory periods (Figure 4) 
mediated by VGSCs (Figures 5 and 6). Our data grant a 



dogma that short-time pulse initiates a single spike at 
axonal hillock [1,4,6], due to low thresholds and high 
density of VGSCs at this segment [2,3,7,9,10,56-67]. It is 
noteworthy that the durations of in vivo fluctuated sig- 
nals are over 50 ms [15] and unitary synaptic events last 
for longer than 20 ms [13,68]. Despite a high density of 
VGSCs at axonal hillock, depolarization pulses above 
these durations make the number of functional VGSCs 
not being high at AIS [26] as well as convert somatic 
threshold and refractory periods to be lower (Figure 4), 
such that physiological signals induce sequential spikes 
being somatic in origin [15]. 

What are physiological impacts for the input-dependent 
plasticity of spike initiation location? This plasticity is re- 
lated to saving the energy during neuronal encoding. The 
neurons step over two barriers [43], spike threshold and 
refractory period [32,44-46], to produce the spikes. The 
energetic barriers to fire the spikes at the soma are 
lower by giving long-time pulses, and these values at the 
axon are lower by short- time ones (Figure 10). There- 
fore, the relocation of spike initiation is used for the 
neurons to save the energy in spike encoding. Moreover, 
the propagation of soma-integrated input signals in 
long-time depolarization toward the axons may inac- 
tivate their VGSCs and be shunted by GABAergic 
receptor-channels at axonal hillock. To prevent the 
decay of these signals in their propagation, the integra- 
tion of analog synaptic signals and the encoding of 
digital spikes are better fulfilled in a single subcellular 
compartment, i.e., neuronal soma. On the other hand, 
fluctuated pulses allow VGSCs to be recovery (Figure 4) 
and even facilitate Navl.6 activation at the axons [69], 
such that the axons serve to an optimal match for the 
frequency between input signals and spikes as well as 
optimal response to high frequent inputs (Figures 1, 2, 3 
and 11). 

We have studied the regulation of spike initiation 
relocation between the soma and axon of cortical pyr- 
amidal neurons, a role of VGSC dynamics in this input- 
dependent plasticity and its physiological impact for the 
neurons to save energy. The studies in these three as- 
pects reveals the natural features of spike initiation and 
strengthens the local presence in the input-dependent 
relocation of spike initiation. 

Methods and materials 

Brain slices 

The study and all experiments were fully approved 
by the Institutional Committee of Animal Care Unit in 
Administration Office of Laboratory Animals Beijing 
China (ID# B10831). The slices from sensory cortices 
(300 [im) were prepared from FVB mice. Mice in post- 
natal day 15-20 were anesthetized by injecting chloral 
hydrate (300 mg/kg) and decapitated with a guillotine. 
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The cortical slices were cut with a Vibratome in the 
modified and oxygenized (95% 0 2 /5% C0 2 ) artificial 
cerebrospinal fluid (mM: 124 NaCl, 3 KC1, 1.2 NaH 2 P0 4 , 
26 NaHC0 3 , 0.5 CaCl 2 , 5 MgS0 4 , 10 dextrose and 5 
HEPES; pH 7.35) at 4°C, and then were held in the normal 
oxygenated ACSF (mM: 126 NaCl, 2.5 KC1, 1.25 
NaH 2 P0 4 , 26 NaHC0 3 , 2.0 CaCl 2 , 2.0 MgS0 4 and 25 
glucose; pH 7.35) 35°C for 1 hour before the experi- 
ments. A slice was transferred to a submersion cham- 
ber (Warner RC-26G) that was perfused with normal 
ACSF for electrophysiological experiments [70]. 

Dual recording 

The soma and axonal bleb [71] of identical pyramidal cells 
in layers IV-V of cerebral cortex were simultaneously 
recorded (MultiClapm-700B, Axon Instrument Inc. CA 
USA) under a fluorescent/DIC microscope (Nikon FN- 
E600). The identification of axonal blebs rather than 
dendritical bleb was based on the diameter and branches 
of processes as well as the polarity of neurons [13,15]. 
Neuronal processes with less branches and fine diam- 
eter are axon. The electrical signals were inputted into 
pClamp-10 (Axon Instrument Inc. USA) with 50 kHz 
sampling rate. 

In whole-cell recording on the soma and axonal bleb, 
the action potentials were induced at these two compart- 
ments, respectively, by long-time steady depolarization 
pulses or cosine waves in various intensities, but the 
stimulus intensities to the soma and axonal bleb of each 
neuron were identical. As the excitability varied among 
the neurons [68], the stimulus intensities for different 
cells were normalized, in which we set the threshold 
current as 1, and then increased pulse intensities by 10% 

gradually, i.e., 1, 1.1, 1.2 2.2 times of threshold 

current. This strategy made the data from different neu- 
rons to be statistically comparable. A judgment for record- 
ing two sites from an identical neuron was based on the 
presence of direct and corresponding electrical signals. 
Transient capacitance was compensated and output band- 
width was 3 kHz. Pipette solution contained (mM) 150 
K-gluconate, 5 NaCl, 0.4 EGTA, 4 Mg-ATP, 0.5 Tris- GTP, 
4 Na-phosphocreatine and 5 HEPES (pH 7.4 adjusted by 
2 M KOH). The osmolarity of pipette solution freshly 
made was 295-305 mOsmol, and the pipette resistance 
was 10-15 Mf! 

It is noteworthy that the axonal blebs of pyramidal 
neurons were formed from the resealing at the end of 
cut axons during slice preparation. Although this prepar- 
ation could be argued as an injured axon, several lines of 
evidence indicated the functional intact of axonal blebs. 
The values of resting membrane potentials and action 
potentials are closely normal (Figures 1, 2, 3, 4 and 5). 
These axons underwent the functional plasticity, i.e., 



upregulation and downregulation in their abilities to fire 
spikes (Figures 8 and 9). Moreover, other studies in 
immunohistochemistry and electrophysiology suggested 
that the functions of axonal blebs were likely normal [7]. 

The intrinsic properties of the somata and axonal 
blebs in our studies included spike thresholds and 
refractory periods measured under a whole-cell current- 
clamp. Spike thresholds were measured by depola- 
rization pulses with the inverse changes in their in- 
tensities and durations. Spike refractory periods were 
measured by injecting two pulses (5% above threshold) 
with various durations into the neurons after each 
spike, in which inter-pulse intervals were adjusted 
[32,33,55,68,72]. In the measurement of spike refrac- 
tory periods versus pulse durations, the increment of 
depolarization duration was inversely associated with 
the reduction of its intensities, which were 5% above 
thresholds. 

Pharmacological studies 

Anemone toxin (ATX; [34] was dissolved freshly in the 
standard ACSF at 5 uM for a final concentration, and 
was given to a specific subcellular compartment by a 
pipette perfusion via the pressure. In the application of 
ATX, fluorescent Alex-488 (10 uM) was also added into 
this solution in order to show the size of areas for per- 
fusing ATX to be restricted to specific subcellular com- 
partment. The pressure was added into the pipettes to 
have ATX/ Alex-488 diffusion within an area of less than 
10 urn diameters, which was seen under a fluorescent 
microscope at 488 nm. 

Single channel recordings 

VGSCs' currents were recorded in a cell-attached config- 
uration by MultiClamp-700B and pClamp-10 at the 
axonal bleb and soma of identical pyramidal neurons. 
Seal resistance was above 10 G£l, and pipette resistance 
was 10-12 MO. Pipette solution contains (mM) 120 
NaCl, 2 MgCl2, 10 HEPES, 30 TEA and 0.1 mibefradil 
[13]. The threshold potentials for VGSC activation were 
measured by adding negative voltage-pulses into the re- 
cording pipettes. 

Modeling 

Computational simulation was achieved in NEURON 
(v7.0), and the following facts were taken into account. 
The dynamics of axonal and somatic VGSCs (Figure 6) 
was introduced into the axon and soma of a neuronal 
model, respectively, to test the initiation of action 
potentials [42,73,74]. Other parameters about VGSCs 
were based on Hodgkin-Huxley kinetics and current 
reports [7,9,10]. 
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The equations for the Na current were based on the 
works of Hodgkin and Huxley. Na current was calcu- 
lated by Ohmic relation: 

l na =g na m 3 h(V m -E na ) 

m and h were calculated by first order kinetic equations: 

, _ (m in f-m) 



h 



(hinf-h) 



Th 



m in f and h in f represented the steady-state value of m 
and h. r m and r h were time constant. m in f and r m were 
calculated by following equations: 



m 



inf 



1 



In these equations, a and /? were the functions of local 
membrane potential (Vm), which were calculated by: 



a(V m ) 



A(y m -Vth) 

l- e -{V m -V th )/K 

-A(V m -V th ) 

l- e -{v m -Vth)h< 



All constants in these equations were shown in Table 
one. In terms of /z, was described analogously to r m , 
but h in f was given directly by: 



h 



inf 



1 _|_ e (V m -V th )/,< 



The constants of this equation also are shown in Table 1. 

The distributions of VGSCs in these two compart- 
ments were based on the references [2,3,7,9,10]. VGSCs' 
reversal potential was set at 50 mV. For potassium 

Table 1 Parameters for the equations describing VGSCs 



Channels 


Variable 


Function 


A(/ms) 


V th (mV) 


K(mV) 


VGSC on soma 


m 


a 


0.182 


-28.2 


9 






P 


0.124 


-28.2 


9 




h 


a 


0.0091 


-50 


5 






P 


0.024 


-75 


5 






hinf 




-55 


6.2 


VGSC on axon 


m 


a 


0.182 


-35 


9 






P 


0.124 


-35 


9 




h 


a 


0.0091 


-50 


5 






P 


0.024 


-75 


5 






hinf 




-59 


6.2 



channels, high-voltage-activated K + channels (Kv) and 
slowly-activated/non-inactivated M-type K + channels 
(Km) were added into this neuronal model. To have 
the initiation of sequential spikes, the densities of Kv 
and Km were 50 and 100 pS/cm 2 on these two com- 
partments, respectively. The reversal potential for K + 
channels was set to -77 mV. In addition, cylinder axon 
was calculated based on 1.6 um in diameter and 70 um 
in length, as well as the soma was 30 um in diameter. 
Their passive electrical properties include C m = 1 uF/ 
cm 2 , R m = 15000 a/cm 2 and R { = 100 O/cm. Resting 
membrane potentials were -71 mV. The simulations 
were run under 37°C. The time step was 0.025 ms. 

Data analyses 

The data were analyzed if the soma and axon had resting 
membrane potentials negatively more than -63 mV and ac- 
tion potentials above 90 mV. Criteria for the acceptation of 
each experiment also included less than 5% changes in rest- 
ing membrane potential, spike magnitude, input and seal 
resistance during each experiment. The data for VGSCs' re- 
cording was taken into account if seal resistance reached to 
10 G£X The values of spike threshold, refractory period and 
VGSC currents are presented as mean ± SE. The statistical 
analyses and comparisons of the results among the different 
subcellular compartments are done by paired t-test. 

Additional files 



Additional file 1: Figure SI. The same inputs to different sites induce 
different outputs, while the propagation of spikes is faithful. The insets 
show simultaneous recordings on soma and axon bleb. The red curves 
are recordings from soma while the blue curves are from axon. A) The 
same long-time step pulse is injected to axon (left panel) and to soma 
(right panel) separately. The outputs are different, but keep consistent 
between the two recording sites. More spikes are induced when long-time 
step pulse is injected to soma (right panel). B) The same fluctuated signal is 
injected to axon (left panel) and to soma (right panel) separately. The outputs 
are different, but also keep consistent between the two recording sites. In this 
case, more spikes are induced when fluctuated signal is injected to axon. 

Additional file 2: Figure S2. Current-voltage relationships of axon and 
soma. A) Membrane potentials of axon (blue) and soma (red) under grade 
subthreshold pulses to these two sites respectively. Schematic for step 
currents is at the bottom of the panel. Dotted lines indicate the points which 
are chose to calculate current-voltage correlation. B) Correlations between 
input currents and membrane potentials of axon (filled circles and blue line) 
and soma (open circles and red line, n=9). There's no obvious difference 
between two lines. 

Additional file 3: Figure S3. ATX does not significantly influence spike 
thresholds and refractory periods at the soma of cortical pyramidal neurons. 
The spike thresholds and refractory periods are measured dynamically by 
changing the patterns of input signals at the soma. A) shows threshold 
stimuli vs. depolarization time at the soma. Compared with the control (red 
symbols), ATX does not change somatic spike thresholds significantly (green 
symbols; n=1 2). B) illustrates spike refractory periods vs. depolarization time 
at the soma. Compared with the control, A~D< does not change the 
refractory periods of somatic spikes (green symbols, n=1 1). 
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